THE GROUND STATE OF RELATIVISTIC IONS IN THE LIMIT 
OF HIGH MAGNETIC FIELDS 



D. H. JAKUBASSA-AMUNDSEN 



Abstract. We consider the pseudorelativistic no-pair Brown-Ravenhall op- 
erator for the description of relativistic one-electron ions in a homogeneous 
magnetic field B. It is shown for central charge Z < 87 that their ground 
state energy decreases according to \/B as B — ♦ oo, in contrast to the non- 
relativistic behaviour. 



1. Introduction 

A relativistic atomic electron of mass to in a magnetic field B = V x A 
resulting from a vector potential A is described by the Dirac operator H \W\ . 

H = Da + V, Da = apA + pm, (1.1) 

where = P ^ eA, a., (3 are Dirac matrices and V — is the Coulomb field 
generated by a point nucleus of charge Z fixed at the origin. The coordinate and 
momentum of the electron are denoted, respectively, by x and p (with cc = |x| = 
^Jx\+ x\+ x\ ), and the field strength is 7 = Ze^ . Relativistic units {h — c—1) 
are used, with « 1/137.04 being the fine structure constant. 

A widely used pseudorelativistic operator which accounts for the spin degrees 
of freedom is the Brown-Ravenhall operator h^^. It can be obtained from a pro- 
jection of H onto the positive spectral subspace of the electron at y = [3] (see 
also [7] for its mathematical analysis). 

The Brown-Ravenhall operator in a magnetic field [H [T31 [TO] is given by 

h^"^ = EA + V1 + V2 (1.2) 

X tiA + TO X tjA + ^71 

where Ea — \Da\ is the kinetic energy operator. 



Ea = \/ (cp/i)^ + "m?- = \J p\ - ea-'B + m? > to, Ae = 

(1.3) 

Here, er = (cri, cr2, (73) is the vector of Pauli spin matrices, h^^ acts in the Hilbert 
space 1/2 (R^) ® and its form domain is i?i/2(M'^) (g) where i?i/2 denotes a 
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Sobolev space. For B bounded or B e L2(M'^) it has been shown [TnHH] that the 
potential Vi + V2 is relatively form bounded with respect to Ea, 

\{ijAVi + V2)ij)\ < ^[i:,EA^) +7cs|l^f (1.4) 
7c 

for V e -ffi/2(K^) ® C^, where 7c = f, cb = f Ve||B||oo, respectively 7^ = 
;| — e, Cb — c(e) ||B|i2 with e > arbitrarily small and some constant c(e) for B 
bounded respectively B S i^(R^). For subcritical potential strength, 7 < 7c, the 
form bound is smaller than one such that h^^ is bounded from below, allowing 
for its extension to a self-adjoint operator. Most recently [13] it was shown that 
for locally bounded A the critical potential strength for the semiboundedness of 
h^^ could be increased to 7c = 2/(^ + |). However, no exphcit B-dependence of 
the lower bound of h^^ was provided. 

In the following we take B — Be^i to be a constant magnetic field along the 
es-axis, generated by 

A(x) - I {~X2,xuQ) (1.5) 

which obeys V • A = 0. 

A matter of interest is the behaviour of the ground state energy when B — + cxi, 
which until now was only rigorously studied for systems with small nuclear charge. 
In this nonrelativistic limit the Dirac operator H turns into the Pauli operator 
Hp = {(TPa)^ + V. For this operator, using a scaling property by which the pair 
{B, ^) changes into {Bq, ^) where Bq is a fixed field and A ^ as S 00, 
the ground-state asymptotics could be investigated by means of a perturbative 
treatment of the electric potential It was derived that, as conjectured earlier 
[TT] . the ground-state energy decreases according to (Ini?)^ as i? ^ 00, the error 
being of the order of In S • ln(ln B) . Recently the Dirac operator itself was studied 
for finite magnetic fields and small nuclear charge (including the simultaneous 
limits _B — > 00 and Z — > 0) in the context of the lowest bound state's diving 
into the negative continuum for sufficiently large -B [5]. A similar scaling as in [2] 
was introduced, and an upper and lower bound for the ground-state energy was 
provided which behave like B^ . 

As concerns the Brown-Ravenhall operator h^^ its ground state was esti- 
mated with the help of a variational wavefunction, similar to the one used in the 
Schrodinger case [TB], with the result that the variationally determined ground- 
state energy decreases asymptotically like B^ jT^]. Recalling that the variational 
energy is an upper bound while the form boundedness (|1.4p provides a lower bound, 
the ground-state energy of h^^ is thus, like that of the Dirac operator, sandwiched 
between bounds which decrease like B^ for B ^ 00. 

The aim of the present work is to prove a stronger statement. 

Theorem 1. Let h^^ — Ea + Vi + V2 be the Brown-Ravenhall operator for an 
electron in a high central Coulomb field of strength 0.1 < 7 < ;| (^^ < 87) and in 
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a homogeneous magnetic field B. Then the ground-state energy behaves like 

Eg c^/B [B oo) (1.6) 

where c > is some constant. 

2. Scaling property 

Let us apply the scaling introduced by [2 to the Brown-Ravenhall operator, 

and set B =: (IqBq with fiQ > and Bq some fixed constant field. Then fiQ oo 

1 

as B oo. Define Xk '■— fJ.§Xk (k — 1, 2, 3) such that 

Pk = -i d^^ = -i fj.^ di^ = Pk 

A(x) = ^(_^-^i2,/,-^ii,0) =:4Ao(i). (2.1) 
Then with m/ ^JJlq —: ifi, 

/i^^ = 4 {Ea + Vi + V2) =: 4 /i^^- (2.2) 



where 



Ea := Jp\-e(TBo +frT?, Pa = P - eAo(x) 



V, := -,A,-A,. v.. := -,A, ^ - ^ ,2.3) 



+ TO i i?^ 

and Ag, follows from p.3p where EA,m is replaced by EA,m. If we introduce 

the scaled ground-state energy iJg := ^ Eg then the corresponding eigenvalue 
equation turns into 

h^^ ^g = Eg i,g. (2.4) 

In order to prove Theorem [T] we thus have to establish two items, 

(i) the existence of a ground state of h^^ for small to > and fixed Bq > 0, 

(ii) the convergence of Eg{fh) as to ^ to the ground-state energy i?g(0) of 
h^^ at TO = 0. 

As a consequence of (ii), 

lim Eg/ 4 - lim Eg{m) - ^(0), (2.5) 

implying Eg ^ Eg{0) — —cB^ as i? ^ oo where c := — i?g(0)/BQ . 
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3. Existence of the ground state 

Lemma 1. Let h^^ be the Brown- Ravenhall operator for an electron of mass m 
in a constant magnetic field Bq. Then for 0.1 < 7 < there exists rUc > such 
that h^^ has a discrete ground state below m for < m < nic- 

Equivalently, h^^ — m has a discrete ground state below for m g 

[0, mc). 

Proof. 

First we show that the spectrum of h^^ is discrete below to, if nonempty. From 
Ea > to we have a{EA) C [to, 00) and hence aessiEA) C [to, 00). Actually 
TO e a ess {E a)- This follows from m G cress (-Da) (when to > 0) for a constant mag- 
netic field jlHl p. 202]. The spectrum is invariant if the unitary Foldy-Wouthuysen 
transformation Uq is applied to Da, which gives UoDaUq^ = (3Ea [3 and proves 
the statement. Since the magnetic field does not restrict the electronic motion par- 
allel to the field, g^Ea) is continuous and hence aess{EA) = [m,co). From ()3.7p 
below it follows that this holds also for to = 0. 

For TO > it was proven in [TU] that aess{h^^) = <^ess{EA)- The respective 
proof, using the compactness of the difference of the resolvents of h^^ and Ea, has 
to be extended to the case to = 0. This is straightforward because the compactness 
property relies basically on the invertibility of the operators Ea+H and h^^ +fj, for 
a suitable n > as well as on the relative form boundedness (|1.4|) of the potential 
Vi + V2 with respect to Ea with form bound smaller than 1. For to = 0, (|1.4p 
remains valid and thus the semiboundedness of h^^. An additional ingredient of 
the proof is the compactness of the operator K :— xo {Ea + where xo is a 
bounded nonnegative function in coordinate space with xo ^ as a; — > cxd. Since 
an error occurred in |10j . we give that proof anew for to > 0. 

Let S"^ := (p — eA)^ -I- to^ be the free Schrodinger operator (increased by 
TO^) in a magnetic field and decompose 

K = [xo {Sa + t^)-'] {Sa + ti) {Ea + A*)"'- (3.1) 

The compactness of the first factor can be shown with the help of the diamagnetic 
inequality. One has the pointwise inequality [I] for any A e L2.ioc{^'^) and V' G 
L2(K^)(8C2, 

l(^V^)WI ^ {l^\m^) (3.2) 



where Ep =^ y'p^ + rri^ . Upon multiplication with xo it follows that |xo('5'a + 
A*)~^V')(x)| < {xo{Ep + \tp\){x). The operator Xo{Ep + is compact as 
a product of two bounded functions of x respective p tending to zero at infinity. 
Therefore Xo{Sa + is compact, too [H Thm 2.2]. 

Concerning the boundedness of the remaining factor in (jS.ip we have for 
|B| < Bo and V £ -H'i/2( 



WSA^Pf = (^, (i^i + e^B) ^) < WEA^Pf + eBoM^ (3.3) 
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whence ||5aV'II < \\EAtp\\ + VeBo With ip := {Ea + Ai)"V it follows that 

1 1 



Ea+ n Ea + ^J. Ea + n Ea + 

< c\M (3.4) 

with some constant c. 

Having thus established that (Jess{h^^) — <^ess{EA) for to > and 7 < ;|, 
we have (Jess{h^^) — [m,oo) under the same conditions. 

Next we show that a{h^^) n {-00, m) = ad{h^^) ^ 0. 

(a) Case m — Q 

In I2J it was shown numerically for relativistic atoms {Z > 20; using the 
scaling and = 1) that there exists a trial function ■0t G i2(K"^) C^, 

V,,(x) = Nt e-^^^^e'l^ g-z.„V-l+i/(ei3o) Q (3,5) 



with Q = xf + X2 and Nt a normalization constant, and an effective charge 
Zef / > such that 

{i}}t,h^^ i^t) < form = 0. (3.6) 

Thus ad{h^^) ^ which assures the existence of a ground state for a given nuclear 
charge Zq. A trial function for Zq satisfying (|3.6p then also satisfies this inequality 
for all Z > Zq. This follows from the fact that the negative part —{Vi + V2) of 
f^BR ^^^^ symmetric operators and hence are < as is V) increases 

linearly with 7 (i.e. with Z) while the positive part Ea is independent of Z. We 
remark that a variational solution to (|3.6p can also be found for Z < 20, however 
the convergence of the respective integrals in {ipt, h^^'ipt) gets increasingly poor 
when Z becomes smaller. 

(b) Case m > 

We consider the operator h^^ = h^^ — m where the subtraction of the rest 
energy m has the advantage that a^ssih^^) = [0, cxd) implying that the ground- 



state energy is below zero for all to. We show that {ipt, h_ ipt) with the fixed ipt 
from p.Sp is monotonically decreasing with m for m < rric when TOc is sufficiently 
small. For the kinetic part we have, using that Ea'i^i — + m? ipt (see also 
section 4) 

{iPt,{EA-m)i;t) = i$tA\/k^ + m^ - m) 0t), (3.7) 
Mk) = 




where (pt G i2(M) is (up to a constant) the Fourier transform of exp(— Ze//A/a;| + I/cBq), 
Ki a modified Bessel function and > a normalization constant [12], [SI p. 482]. 
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With the mean value theorem apphed at to = we get 

{^t,{^jB+m-^ -m) 4>t) = {4>u\k\4>t) - m-Dt{0 (3.8) 
with ^ G [0, to] and 



DtiO := 2iV^ / dk 5^-2 , — (3.9) 







Since the second factor in (|3.9p is nonnegative, bounded by 1 and decreasing with 
^, the integrand is positive except for /c = leading to < Dt{mi) < Dt{^) < 1 
for all ^ and to below some fixed value toi. Consequently, (j3.8p is decreasing with 
TO (linearly for to — > 0). 

Fourier transforming the expectation value of the potential part of h^^, we 

have 

poo poo 

{A,{Vi+V2)i't) = -7C0 / dk / dk' Mk)Vo{k~k') fkk'im) Mk') (3.10) 



A.(TO) A.(fc)A.(fc') + A^ik) A^ik') 

where cq is a constant, EA{k) — + ^ti^, A£;(fc) = ( ^2]t}(k)^ ) ^ ' ^^"^ ^(9) is the 
Fourier transformed expectation value of - with the transverse function (^^'^Boq /4 
(for details see section 4, in particular ()4.6p ). 

One can split F{m) :— {tpt, {Vi + V2) V't) into F+(to) + (to) where F+{m) 
results from kk' > in the integrand of p.lOp while F_ (to) is assigned to kk' < 0. 
The function fkk' is positive and analytic for to, > if (fc, k') e M.'^\Sf with 
Sf (0,M) U (M,0). Concerning F+(to), i.e. kk' > 0, one gets by elementary 

computation %^\m=o = and ^-^^\m=o ^ (except for k = k'). Thus /fet' 
has (almost everywhere) a local extremum in to = 0. Since the remaining factors 
in the integrand of (j3.10p are independent of to and positive (see Lemma [5] for 
Vo), F+(to) has also a local extremum in to = 0. 

Thus there is an TOc such that the variation of .^^(to) in [0,TOc) is consid- 
erably weaker than the linear dependence p.Sp of the kinetic term. The other 
contribution, F_(to), is monotonically decreasing like the kinetic term. In fact, if 

kk! < 0, fkk' is monotonically increasing with to (as 1 - (E^{k)+m)'{EA(k')+ni) is 
increasing with to). As a result, {ipt, [h^^ — m) tpt) is decreasing with to in [0, TOc). 
This feature is confirmed numerically, with a quite large TOc. 

Consequently, for a given nuclear charge Zq, once {ipt, ipt) < is es- 
tablished at TO = it follows that {tpt, tpt) < for to G [0,TOc). On the 
other hand, as discussed in the context of to = (but remains true for to > 0), 
tpt) < for all charges greater than Zq. Since Zq > 10 is arbitrary, this 
guarantees the existence of a ground state of h^^ below for all to < nic- □ 

Alternatively, the existence of a discrete ground state (and, in addition, infin- 
itely many bound states) of h^^ below to (for to > 0) may be based on a theorem 
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by Matte and Stockmeyer [15J . However, its applicability requires the fulfilment of 
certain nontrivial conditions on the Weyl sequences for the essential spectrum of 
the free Dirac operator Da- 

In the remaining part of the present work we prove the convergence of the 
sequence of eigenvalues Eg{fh) for to ^ to EgiO). We start by restricting the 
ground-state function to the lowest Landau level and show continuity of the ex- 
pectation value of h^^ at to = for a certain class of functions. Then we allow 
for the presence of higher Landau states. From the continuity property of h^^ the 
convergence of Eg{m) to Eg{0) is deduced. 



4. Reduction to a one-dimensional problem 

We want to gain information on the ground-state wavefunction of h^^ for 
arbitrary mass to > 0. To this aim we reduce the three-dimensional problem 
to a one-dimensional one by invoking the eigenfunctions of the Pauli operator. 
The resulting eigenvalue equation is then transformed to an integral equation in 
Fourier space from which the basic properties of the ground-state function can be 
extracted. 

The eigenfunctions ipnidsixi, X2) = fnidixi, X2) Xs of the two-dimensional 
Pauh operator obey [IHl p. 196] 

leNo, n e No, s = ±1, = [ ±1 n>l ' (^-^^ 

where the spin functions are x+i = (o)j X-i = (i)- These eigenfunctions form a 
complete set of orthonormal functions (see e.g. [21 §111], El) such that an eigen- 
state V' S L2(IR'^) (Xi of h^^ can be expanded in terms of these functions. When 
B strongly dominates the Coulomb field the ground state of h^^ is approximately 
characterized by the ground state of the Pauli operator which is determined by the 
quantum numbers n — and s = +1 (where the rhs of ()4.ip attains its minimum 
zero) . Indeed it was shown for the Dirac operator in a homogeneous magnetic field 
[5] that the n = approximation becomes exact when the nuclear charge tends to 
zero or equivalently, when i? — > 00. 

With the restriction to n = and s = +1 the (normalized) ground state ip 
of h^^ can be written as 

00 

^„=o(x) = J2aiNi e'-^-B-'l^ e^^^ cfiixz) X+i, 
1=0 



(eBo 
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where a/ is an expansion coefficient, g = y/x\+x^, 4'i{x3) ^ normalized function 
yet to be determined, and we fiave used the exphcit form of i^oz ilSj p. 196], [HI the 
negative sign of d in that work relates to a positive coupling in ()4.1|) . {pi + eAi)]. 
When n — and the spin is fixed in ip one can show by using crpji^crpA = 
Pa-^Pa + "icriPA X ■^Pa), that not only Ea and Vi, but also V2 is diagonal in I. 
For s = +1 we have from (|4.ip for the kinetic part {tpoi, Ea ^poi) — \/p^+m^ —■ 
Ea{p3) and it is straightforward to verify for the potential part 

Vm{x3,l) ■■= ix+i foh-- {Vi + V2) ipoi X+i) 
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Ae{P3 



i+i 



2ni 



dg g' 



P3 



21+1 g-eSo£)V2 



P3 



(4.3) 



AEiP3) 



^ + x\ EAiP3) + m y/g-^ + xj EA{p3)+m 

where the functional dependence of Ae{p3) on Ea{p3) is given by (|1.3p . The small- 
est value of {ip,h^^tp) is given hy I — 0. In fact, £^^1(^3) is independent of I. 
Moreover, the weight factor of {g'^ + Xg)"^ in the integral in (|4.3p is peaked at 
gmax = \/ (2/ + l)/ei?o, its normalized value at g^ax slightly decreasing with I. 
Thus Ymi^X'i, I) attains its maximum at I = where gmax is smallest. For the consid- 
eration of the ground state we can therefore restrict ourselves to ^ = 0. Thus, drop- 
ping the sum in (|4.2p and the index on (jj, we have ■ipn=ai^) ^ ^'00(2:1, X2)(f>{x3)(^) . 
Disregarding the coupling to higher Landau states, the eigenvalue equation for 



7 Vmix3,0) 



(t>ix3) = Em- 0(X3). 



(4.4) 



where Em- is the ground-state energy of h^^ — m under the restriction n = 
(and s ~ +1). 

When Fourier transforming (j) in (|4.4p and projecting with an eigenstate to 
Pa, exp(ipa;3), the momentum operators in Vm turn into functions of k and p, 
respectively. As a result we obtain an eigenvalue equation for the momentum- 
space function (/), 



\\/ P^ + m"^ — m) (j){p) j= / dk Ae{p) 



Vo{p-k) + 



Ea{p) + m 



Vo(p-k) 



EAik)+m 



AE{k) (l>{k) = Em- <t>{p) 



(4.5) 

for p G M, where we have introduced the Fourier transformed expectation value 
^0(2^3) of i p.419]. 



dx3 e / gdg e' 
'0 
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2eBQ 



edBe-^'B<^^''^ K^{\q\g). (4.6) 

Lemma 2. The momentum- space potential Vq from \4.6^ obeys 

{-ciln|g|, 
(4.7) 

with ci > and C2 > 0. Va{q) >0 is monotonically decreasing with \q\. 

Proof. For q the modified Bessel function diverges logarithmically, iiTodglfp) ^ 
— ln(|(7|p) = — In \q\ — In g, and in this limit the integral over g is convergent. For 
large [gj we make the substitution \q\g = z and obtain 

where we estimated the exponential by unity and used [21 p. 684]. Since ifo is 
monotonically decreasing and the integrand > 0, so is Vo((7). □ 

With the help of Lemma 2 some properties of (t>{k) can be derived. We have 

Lemma 3. For m > the restricted ground-state Junction in momentum space, 
<j), has the following properties 

(a) 0(fc) is uniformly bounded for /c G M. 

c 

(b) 10(^)1 < 1 — r, \k\ — > oo, where c> is some constant. 

\k\ 



Proof. For the proof of boundedness we use that \/p^ + rn^ — m+ \Em- \ > \Em-\ 
and estimate from (|4.5p . using that Ae < 1 and k/{EA{k) + m) < 1 (for m > 0), 

|0(p)| < ^ — ^ dfcl/o(p-fc) |0(fc)| 



< 



/2^ 
27 1 



(^J^ dq {Vo{q))y II0II, (4.9) 



/2^ |S„ 

where the Schwarz inequality was applied. Due to Lemma 2 the integral is bounded, 
and \\(j>\\ — 1, such that \(l){p)\ < — j- where the constant c does neither depend 

on (j) nor on m. Recall that Em- < exists for all m > since the trial function 
is an n = (spin-up) state. 

It is now easy to derive the behaviour of 4>{p) at \p\ — > 00. Without restriction 
we can assume \p\ > m such that \/p'^ -\- rrfi — rn + |i?m-| > . ^ > . 

YP^+m^ +m v2+l 

For the estimate of the integral in (|4.5p we can use ()4.9p . Then, with the same 
constant as above, 

\4>{P)\ < - c for|p|^c3o. (4.10) 



10 



D. H. JAKUBASSA-AMUNDSEN 



□ 

5. Continuity property of h^^ in m = 

Let 

Mpos W £ -^2(1^) : = 1, uniformly bounded near 0} (5-1) 

be the set of square integrable functions for which there exists C>0, 0<(5<1, 
such that |<y3(A:)| < C for k e [—(5,(5]. Then for any eq > there exists < (5o < (5 
with 

2 (sup \^\f-26o < AC^So < e^. (5.2) 

[-s,s\ 

We can choose 5^ = min{;j|^, 5} for aU Lp e Afpos- 

Lemma 4. Let h^^{m) be the Brown- Ravenhall operator for a particle with mass 
m and let iP{k) = Nq 6-"^°^^ ^'^(t){x3) x+i with (j) e Mpos- Then {ip, h^"{m) t/j) is 
continuous in m = 0, i.e. for any e > 3 mo > : 

KV', (/i^^(m) - /i^^(O)) V)l < e for all m < tuq. (5.3) 

As a consequence, the expectation vahie of h^^{m) is also continuous in 0. 

Proof. Denoting the dependence on m by an additional subscript we have for the 
kinetic part of h^^{m) from (|1.3p . 

9 

77? 

\{rP,{EA,n-EA^o)i^)\ < -J^^W < (5-4) 

which shows the continuity in m = 0. Note that (|5.4p holds for all ifj E L2{M.^)<E)C^ . 
The potential part can be decomposed in the following way, 

Vl,m - Vi,o = -J As.m- {A-E.m - Ae.o) " 7 (^B,m " ^B.o) " ^E.O (5.5) 

X ' X ' 

with a similar formula for V2. The expectation value of the second summand is 
the conjugate complex one of the first factor (except for the replacement of AE,m 
by Ae.o) and need not be considered separately. Transforming into Fourier space 
we have from (|4.6p 

\{lp, {AE,m - Ae,o) - AE,,n V")! = \{<t>, (AE.niip^) " ^£,0(^3)) Vo AE,m{p3) 
X 

/oo />oo 
dk\4>{k)\ dk' K{k,k')\4){k')\ (5.6) 
-00 J —00 

with the kernel 

K{k, k') \AE.Ak) - AE^m Vo{k ~ k') AE,^{k'). (5.7) 

v 27r 

We estimate (j5.6p further by applying the Schwarz inequality, 

KV', {AE.m - Ae.^)\Ae,^^)\ <([ dk\${k)\^Iik)\ ' ( I dk'\4>{k')\^J{k')\ ' 
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/(/c) := / dk' K{k,k'), J(k'):^ / dkK{k,k'). (5.8) 



We will show that each of the two factors can be bounded by an arbitrarily small 
£2^ provided m is sufficiently small. We have, substituting ^ ~ k — k' and using 
that AE,m{k') is bounded by 1, 

I{k) < -^\AE,,n(k)-AE.o{k)\ I diM^), (5.9) 

V 2,11 J-oo 



where by Lemma [2] the integral is equal to some finite constant c/. We make use 
of the fact that AE,m{k) is continuous at m = if fc 7^ 0, choose eo > and take 
(5, 5q from the definition (|5.ip of Mpos- Then for |fc| > 60 there is moo > such 
that \AE,m{k) — AE,o{k)\ < eo for all m < moo- We decompose the integration 
interval according to (— cx),oo) = (— 00, — (5o) U [— (5o,(So] U (Sq^oo) and obtain 

dk \m\^ m < ^\eo( f '\ r]dk \m\' 

, y/Zir yj-00 JSo J 

So ] 

dk |0(fc)p \AE.mik) - AE.oik)\ > ■ (5.10) 

-So J 

We have \AE,m{k) — AEfi{k)\ < 2, and since (j) G Mpos, by ()5.2p the last term in 
the curly brackets is estimated by 2 ( sup |(/)| )^ • 2So < eo- The other two terms 

[-<5o,<5o] 

can be estimated by eo H^IP = eo. This leads to 

2cr 

dk |(/>(fc)p I{k) < -jL eo for m < m^oo- (5.11) 
! V 2vr 

For J{k') we treat the case \k\ > Si with 5i < 1 by estimating \AE.m{k) — 
AE.oik)\ < eo for, say, m < moi. Then we obtain, substituting ^ k ^ k' for k, 

J{k') < -L AB,„(fc')|eo ' "^^^"^^^ 

Si~k' ^ 

d^ \AE,mi^ + k') - AEflii + fc')l MO \ ■ (5.12) 

-Si-k' J 

The first term in the curly brackets is estimated by c/eo. As concerns the second 
term we profit from the fact that Vb(C) is symmetric and monotonically decreas- 
ing with |^|. Thus, estimating the difference between Ae^tu and Ae^ by 2, the 
remaining integral over Vo(^) has its maximum value for k' = 0. With (|4.7p we 
have 

/(Si— fc' <.(5i 
d^MO < / d^VoiO 
-Si-k' J-Si 

pSi 

< -2ci / \n£_d£_ = 2ci (5i(l-ln(5i). (5.13) 
Jo 
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The rhs tends to zero as i5i — > and thus can be made smaller than eo for suffi- 
ciently small ^1. Then, with AE,m ^ 1, 

dk' |0(fc')P J{k') < -^{c/eo + 2eo} for m < moi. (5.14) 

The estimate of the second potential part, V2,m — V2.0 proceeds in the same 
way. One only has to replace throughout AE.m(k) by AE.m{k) := AE,m(k) ■ 
3 which is also bounded by 1 and continuous at to = for fc ^ 0. This 
proves that the potential part is bounded by e for to < min{TOoo, toqi} —: toq. □ 

We continue by showing that the state iprm defining the ground state V'o,™ 
of both h^^{m) and h^^{m) (under the restriction n = 0, s = +1; the subscript 
TO is added for clarity), is a member of Mpos for to sufficiently small. For m = 
one has 0o S Mpos because from (|4.9|) . |(?!'o(p)| < for all p G M (where Eq is 
the ground-state energy of the restricted h^^{0)), and we may choose any bound 
C > Therefore, from LemmaH](with possibly a slightly smaller toq), 

\{^Pofi,h^'^{m)^Pofi)-{i^ofi,h'''^iO)ijofi)\ < e for to < toq. (5.15) 

As a consequence, (V'0,0, h^^{m) ■00,0 ) < Eq + e. 

On the other hand, E^- < {'4'o,o,h^^'{m)iljo^o) since the expectation value 
of (the restricted) h^^{m) taken with an arbitrary function leads to an upper 
bound of the ground-state energy Em- ■ Combining these two inequalities we get 

Em- < Eq + e for to < toq. (5.16) 

Since _Eo < there exists 62 < such that Eq < S2 < 0. Let e be so small that 
Eq + e < 82- Then, from (|5.16p . Em- < S2 for all to < toq. This leads to an 
TO-independent bound of (pm from (|4.9p , 

c c 

\(t)mip)\ < tt; r < TTT for TO < Too. (5.17) 

\Em-\ M 

With the choice C := > -i^-r we have found a universal bound on chm for all 

\02\ I-C/Ol 

TO < Too- As a consequence. Lemma U] holds for all cigcnstates ipa,m with to < toq. 

Let us now generalize the ansatz for the ground-state function tp of h^^ by 
including the higher Landau states (n > and s = —1) in the expansion (|4.2p . 
This allows for the coupling of the potential Vi + V2 to n > 1 states as well as for 
spin- flip induced by V2. However, for high magnetic fields the ground state is not 
much affected. For the Dirac operator it was shown numerically [H] that for e.g. 
Z = 68 and B = 2.3 x lO^^ G (where m'^e^ = 1.1 x lO^^ q jg ^^e particular field 
for which Coulombic and magnetic forces on the electron are equally important) 
the effect of considering the n > contributions is about 3 percent. 

With tpn=o being the ground state of h^^ with the coupling to higher Landau 
states switched off, we now add a spin- flip term ip-i and a remainder ^pr, 

V'(x) = aoV'n=o(x) -I- ai'(/'-i(x) -I- 02 '(/'r(x). (5.18) 
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ip-i is characterized by 71 = 0, spin s — —1 and I — 1 since V2 changes I by one unit 
simultaneously with the spin- flip, tpr is composed of Landau states with n > 1. All 
functions in ()5.18p are normalized and mutually orthogonal. The weight factors 
obey |aoP + |aip + |a2p = 1, guaranteeing that ip is normalized too. 

Since Ea, acting on a Landau state with n > 1, is according to (|4.ip given 
by Ea — ^2neBo + eBo{l — s) + p§ + to^ > ^2eBo + m? it is strictly positive 
at m = 0. The same is true when Ea acts on an n = spin-down state, resulting 
in Ea = y^2eBo + + m? > \/2eBo + m?. Hence, the m-dependent factors in 
V\ and V2, Ae and (TI>a/{Ea + m), are analytic in m = 0. As a consequence, 
{ipr, h^^ipr) and {-ip-i, h^^'ip-i) are continuous in m = (in the sense that ipr, 
respectively tp-i, is kept fixed when performing the limit to ^ 0). 

When proving Lemma [H for ip from (|5.18p it thus remains to show that 
the off-diagonal matrix elements are also bounded by e for to sufficiently small. 
Considering the potential Vi we have, using the decomposition ()5.5|) . 



-- (V'n=0, - Vl,o) V'r) = {\ Ae ,m i'n^O , {A-E .rn - Ae ,o) i^r) 

7 a:2 x2 

+ {\ {Ae.jh ~ AE,o)'4'n=0, AEfii^r)- (5.19) 
X2 x^ 

The rhs of (|5.19p can be estimated from above by 

\\^Ae,„i V'ri^oll • 11^ {AE,in - ^B.o) 
+ {i'n=0,iAE,m- Ae,())-Ae i^n=o)^ ■ {A, --Vlfi^pr)^ (5.20) 

X 7 

with the bounded operator Ae ■— Ae.tu — Ae^- Since all Landau states, and 
thus ipr, have a Gaussian decay (^ g-eSoe /4 jgj^^ (V'n^i.o^r) is bounded. The 
multiplication factor is bounded by e5 according to the proof of Lemma SI The 
two factors in the first term of (j5.20p are bounded by e since {il)n=o, yi,m V'n=o) is 
bounded and AE,m continuous in to = when acting on 

The corresponding estimates for the potential V2 are done in the same way. 
All these estimates also hold when 4'r in (|5.18p is replaced by V'-i or when ipn=Q 
is replaced by Together with (|5.4p this establishes the continuity property of 
fiBR^ respectively h^^, 

KiPAh^^'im) - h^^{0))'iP)\ < £ forTO<TOo (5.21) 

and Too sufficiently small, where ip from (|5.18p is taken as eigenstate to h^^{m') 
for any fixed to' with < to' < toq. 



6. Convergence of the sequence of eigenvalues 
With (|5.2ip at hand it is easy to prove 
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Lemma 5. Let Eg^{m) he the ground-state energy of h^^{m) — h^^im) ~ m, 
and let Eg{0) be that of h^^{0). Then for a given e > there is an toq > such 
that 

\Eg^{m) - Eg{0)\ < e for all m < mn. (6.1) 

Proof. Let us choose tp in (|5.2ip as eigenstate ip„i to h^^{m). Then 

\{ij^,h^'^{m)i:,n) - {'^m,h'''^{0)^^)\ < e for to < mo, (6.2) 

such that (V'm,ft-^^(0)Vm) < Eg^{m)+e. Moreover, one has £'g(0) < (ipm, h^^iO) ipm) 
since ipm wih differ from the ground-state function ipQ of h^^{0). The combination 
of these two inequahties leads to 

Eg^{m) > Eg{0) ~ € for TO < Too. (6.3) 

On the other hand, with t/j := ipo, the step from (I5.15|) to (j5.16|) can be mimicked. 
Then l|5.16p turns into Eg^{m) < Eg{0) + e which gives, combined with (|6.3p . 

\Eg^{m) ^ Eg{0)\ < e for to < toq. (6.4) 

This proves the continuity of Eg^{m) (and of Eg{m) as well) at to = 0. □ 
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